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ABSTRACT 

We investigate a "hypernova" model for gamma-ray bursts (GRBs), i.e., massive C+O star 
model with relativistic jets. In this model, non-thermal precursors can be produced by the 
"first" relativistic shell ejected from the star. Main GRBs are produced behind the "first"-shell 
by the collisions of several relativistic shells. They become visible to distant observers after the 
colliding region becomes optically thin. We examine six selected conditions using relativistic 
hydrodynamical simulations and simple analyses. Interestingly, our simulations show that sub- 
relativistic (v ~ 0.8c) jets from the central engine is sufficient to produce highly-relativistic (T > 
100) shells. We find that the relativistic shells from such a star can reproduce observed GRBs with 
certain conditions. Two conditions are especially important. One is the sufficiently long duration 
of the central engine > 100 sec. The other is the existence of a dense-shell somewhere behind 
the "first" -shell. Under these conditions, both the existence and non-existence of precursors, and 
long delay between precursors and main GRBs can be explained. 

Subject headings: gamma rays: bursts — gamma rays: theory — supernovae: general 



1. Introduction 

Long Gamma-Ray Bursts (hereafter just GRBs) 
are likely produced by massive stars, because they 
appear in star-forming regions (e.g., Vreeswijk et 
al. 2001; Gorosabel et al. 2003); and some of them 
are followed by "bumps" in the afterglows that 
can be most convincingly explained as bright su- 
pernovae (e.g., Reichart 1999; Bloom et al. 2002). 
Most clear evidence was seen in GRB030329 whose 
afterglow turned into a supernova (SN), SN2003dh 
(e.g., Stanek et al. 2003; Hjorth et al. 2003). The 
spectrum of this SN showed that it was an ener- 
getic type Ic supernova, called "hypernova" (HN, 
e.g., Nomoto et al. 2004). Prototype of a HN is 
SN1998bw (Galama et al. 1998). The bumps in 
the afterglows are also consistent with the light 
from HNe. 



1 Department of Astronomy, University of 
Tokyo, Hongo, Bunkyo-ku, Tokyo 113-0033, Japan; 
umcdaOastron.s.u-tokyo. ac.jp, tominaga@astron.s.u- 
tokyo. ac.jp, nomoto® astron. s . u-tokyo .ac.jp. 

2 Department of Earth Science and Astronomy, 
Graduate School of Arts and Science, University of 
Tokyo, Komaba, Tokyo 153-8902, Japan; maeda@esa.cu- 
tokyo.ac.jp 



Since 1998, four HNe (SN1997ef, SN1998bw, 
SN2002ap, SN2003dh) have been identified through 
the light curve and spectral modeling. The pro- 
genitors of these four are C+O stars of 5 to 14M© 
with explosion energies ~ 4 — 50 x 10 51 ergs (e.g., 
Nomoto et al. 2004). Such C+O stars can be 
formed from massive stars with main-sequence 
masses of M ~ 20 — 40M© after losing outer 
envelopes by the strong wind mass-loss or interac- 
tion with binary companion stars. We note that 
the existence of He layers may not affect on the 
observed properties of the SNe. 

GRBs are produced from the shocks where 
highly relativistic shells collide with other shells 
forming internal shocks, or with circumstellar 
matter forming external shocks. Past studies have 
revealed that the break in the afterglow light curve 
strongly suggests that the relativistic flow is in the 
form of jets with opening angle of 5 ~ 10° (Frail 
et al. 2001). Variabilities in the main-GRB favor 
the internal shock model than the external shock 
model, though afterglows are well-explained by 
the external shocks. 

In this Letter we test a hypothesis that GRBs 



originate from massive C+0 stars, and investigate 
whether the time scales and observed properties of 
GRBs, including precursors, can be consistently 
explained. We expect that adding He envelopes 
to the stars does not affect the results much, be- 
cause He stars have radius only few times larger 
than C+O stars. There are some previous works 
that considered the ejection of a relativistic jet 
from a massive star, but they mainly focused on 
the jet formation in the magnetohydrodynamical 
processes (Mizuno et al. 2004; Proga et al. 2003) 
or the jet propagation inside star and relatively 
soon after the breaking out of the star (Aloy et al. 
1999, 2002; Zhang, Woosley, & Heger 2004), and 
have not discussed much about the properties of 
GRBs themselves. 

We discuss six conditions for successful mod- 
els of GRBs by showing some results of numerical 
simulations and simple analyses to lead to general 
conclusions and to clarify remained issues. More 
detailed models and unconsidered effects such as 
rotation, MHD effects, and the effects of accom- 
panied HNe are further explored in the following 
papers. 

2. Hydrodynamical Models 

We have developed a multi-dimensional special 
relativistic hydrodynamic code (Tominaga et al. 
2005). In this Letter we simulate relativistic jets 
from a 14M C+O star, embedded in the circum- 
stellar matter (CSM) with a density structure r~ 2 . 
The stellar radius is about R = 3.5 x 10 10 cm and 
CSM density at r = 10 11 cm is 10 _11 g cm" 3 . The 
jet is initiated at r = 2 x 10 9 cm with opening 
angle of 15°. At this location, the energy is in- 
jected mostly as thermal energy (E t h/E]^ n ~ 40) 
with v = 0.8c and E ~ 5 x 10 51 erg s" 1 and the 
jet duration time, Atj et = 9 sec. In this model the 
total energy of 4.5 x 10 52 erg is injected as a jet. 

Here we show only the results of 2D calcula- 
tions with axisymmetry and with relatively coarse 
meshes, because the purpose of this Letter is to 
make qualitative discussions. The adopted mesh 
size is 3000 x 45 meshes with more meshes in the 
radial direction in Eulerian spherical coordinate. 
This mesh size is sufficient to estimate the maxi- 
mum Lorentz factors, T, within the error of 10%. 
Higher resolution results as well as other parame- 
ter dependencies, will be explored in the following 



papers (Tominaga et al. 2005). 

The injected sub-relativistic jet is accelerated 
up tp T ~ 50 before breaking out of the star by 
the conversion from thermal to kinetic energies. 
We note that the jet injection has to continue at 
least until the shock front breaks out of the star. 
Otherwise the shock is strongly decelerated by the 
outer envelope of the star. The evolution of the 
peak T is shown in Figure la after t = 9 sec when 
the jet injection is terminated. This figure shows 
that the peak T continues to increase to T ~ 160 
until t ~ 200 sec because of the thermal to the ki- 
netic energy conversion, or thermal expansion. In 
this model, 92% of the total energy is contained 
within 9° from the jet-axis. The highest Mach 
number in the jet is about 10 that is achieved near 
the shock front of the jet. Our calculations are 
similar to those in Zhang et al. (2004). They as- 
sumed narrower and higher velocity jet injection, 
smaller E^ot an d E tot than our models, but the re- 
sults are qualitatively consistent. We can adopt a 
lower velocity for the injected jet to obtain highly 
relativistic jet, because our injected jet is thermal 
energy dominant. 

3. Six conditions to satisfy 

Here, we will examine the following six condi- 
tions for a successful GRB model to satisfy, (i) 
Can highly relativistic jets (r > 100) with suf- 
ficiently large energy E > 10 51 erg be produced? 
(ii) Where should GRBs be produced? Is the GRB 
producing region optically thin? (iii) Can the time 
scale of variability, t var i ~ 1 — 20sec be repro- 
duced? (iv) Can the duration, idur ~ 100sec of 
main-GRBs be reproduced? (v) Can the gener- 
ation of precursors be reproduced? (vi) Can the 
long delay time between the precursor and main- 
GRBs, tdci ~ — 200 sec be reproduced? 

(i) The GRB spectra are non-thermal and con- 
sidered to be produced by synchrotron radiation. 
According to the standard fireball model, relativis- 
tic shells with F > 100 may produce such spectra 
(e.g., Piran 2005)7" 

The relativistic shell must have sufficiently 
large energy to explain the observed flux. For 
some gamma-ray bursts the estimated isotropic 
energy is as large as ~ 10 54 erg. However, in the 
jet model the estimated energy can be as low as 
~ 10 51 erg. Therefore, the jet should have energy 
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at least <~ 10 51 erg, though more than ~ 10 52 erg 
might be better considering energy convergent ef- 
ficiency from the kinetic energy to radiation. 

There are mainly two ways to eject relativistic 
materials from massive stars. One is to eject point 
like masses with relatively high-density and small 
cross sections, known as the cannon-ball model 
(Dar & Rujula 2004). The other is narrow con- 
tinuous relativistic jets that last at least until the 
front end of the jets breaks out of the star (Aloy et 
al. 1999; Zhang etal. 2004). The former type ob- 
jects cannot be significantly accelerated after the 
ejection, thus highly relativistic matter with suf- 
ficiently large density have to be produced near 
the central engine (Umeda 2000). Since very little 
is known for the central engine, we can not ex- 
clude such a possibility, but the latter type proba- 
bly has less obstacles because even sub-relativistic 
jets (v ~ 0.8c near the center) are sufficient to 
produce highly relativistic jets (Fig.l). As shown 
in Figure 1, jets from a C+O star can have suffi- 
ciently large T and E. 

(ii) We have shown that jets from a massive 
C+O star could potentially produce GRBs. How- 
ever the particular model in Figure 1 cannot re- 
produce the observed GRBs because it is the single 
smooth jet and cannot produce variability unless 
CSM is quite inhomogeneous. It is usually consid- 
ered that the variability in GRBs is explained by 
the internal shocks produced by the collisions of 
several shells. Such shells might be produced by 
the variable activity of the central engine. Then, 
several relativistic shells with various T may col- 
lide to produce GRBs through the internal shocks. 
Other mechanisms to produce several shocks have 
also been proposed in the literature. For example, 
Aloy ct al. (2002) suggested that a shear-driven 
instability leads to rapid fluctuations that produce 
the internal shocks. Although further study is cer- 
tainly needed, we may say that the C+O stars 
with relativistic jets are quite promising systems 
to explain the observed variability of GRBs. 

We assume collisions of relativistic shells pro- 
duce main GRBs. Since the colliding regions are 
located behind the first relativistic shell, the GRBs 
are not observable if the first shell is optically thick 
in gamma-rays. In Figure la we show the gamma- 
rays photosphere, where the gamma-ray opacity 
is assumed to be 0.03 cm 2 g _1 . This figure shows 
that around t — 1600 sec, the photosphere passes 



through the "first" relativistic shock. Therefore, 
GRBs are visible only after 1600 sec in the local 
time. Is this too long or sufficiently short? The 
duration of the main-bursts is typically larger than 
10 sec in the observer's frame, that means the du- 
ration is ~ 2r 2 x 10 sec > 2 x 10 5 sec in the lo- 
cal time. Therefore, the optical thickness is not a 
problem. 

(iii) The observed time scale of variability is 
typically i var ; ~ 1 — 20sec in the observer's time. 
In our model the distance between the internal 
shocks determines the time scale. Let / be the 
mean distance between the internal shocks in the 
local frame of the star, then the mean interval of 
two GRB peaks from these shells is l/c for ob- 
servers (Fig. 2). Therefore, the observations re- 
quire that I ~ c x t vari ~ (3 — 60) x 10 10 cm. 
Again, whether this is possible or not depends on 
the properties of unknown central engine. Success- 
ful model should explain this length scale. This I 
has the same order to the radius of the progenitor 
star and this fact may be suggestive. 

(iv) The duration time of the GRBs is typi- 
cally tdur ^100 sec. This requires that the re- 
gion where the internal shocks distribute should 
extend to the distance L ~ c x t<jur ~ 3 x 10 12 
cm in the local frame (Fig. 2). Roughly speaking 
this L should be related to the jet duration time 
Atjot by L ~ c x Atjot, however as described be- 
low this L can be larger than c x Aij et by a factor 
of a few, considering the deceleration of relativis- 
tic shells by the collisions with other denser shells. 
Thus we may write this condition as 

Ai jet > idur ~ 100 sec. (1) 

Without knowing the mechanism of the cen- 
tral engine, it is difficult to discuss whether such 
a jet duration time is reasonable or not. How- 
ever, similar situation has been investigated in 
the "collapsar model" (MacFadyen, Woosley & 
Heger 2001). They suggested a possibility to cre- 
ate narrow jets by the unspecified processes, such 
as the Blandford-Znajck (1977) mechanism and 
other MHD precesses (e.g., Blandford & Payne 
1982). 

It is reasonable to assume that Atj Ct is governed 
by the time scale of mass accretion onto the central 
black hole. In the collapsar model, the accretion 
time scale is short, ~ lOscc, for the Type I case, 
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where the outgoing shock fails to be launched from 
the collapsed iron core. On the other hand, it is 
relatively long, ~ 30 — 3000 sec, for the Type II 
case, where the inner layers of the star initially 
move outward but lack adequate momentum to 
eject all the matter. Therefore, the Type II collap- 
sar model at least can satisfy the above constraint 
(!)• 

In the jet-like explosion model of hypernovae 
(Maeda & Nomoto 2003), most materials are 
ejected along the broad "jets" , while large amount 
of matter may fall-back for other directions. For 
the equatorial direction, the fallback time scale 
may be as long as a few hundred seconds as in the 
Type II collapsar model. 

(v) The precursors are important to constrain 
the models of GRBs. However their properties are 
still uncertain mainly because different authors de- 
fine them differently We follow the definition by 
Lazzati (2005). Most precursors have non-thermal 
spectrum and they have equal softness or softer 
spectrum than the main-bursts (Lazzati 2005). 
They contain a small fraction (0.1 — 1%) of the to- 
tal event counts. The delay time between the pre- 
cursors and main-GRBs, idur, is typically the or- 
der of hundred seconds. This Lazzati's precursors 
may be different from the X-ray precursors, that 
are by definition not necessarily separated events 
from the main GRBs. If the softness varies during 
the main GRBs, earlier or later part of the bursts 
may be identified as X-ray precursors or postcur- 
sors in some definitions (e.g., Nakamura 1999, for 
a model to explain them). 

The existence of the Lazzati's precursors is not 
well established, bu if they are confirmed, long 
delays and the non-thermal spectra are difficult 
to reconcile with conventional GRB models. For 
example, if both events are produced by the same 
engine, it implies that after the engine forms a 
precursor, the activity is silent for 100 seconds and 
then main-activity restarts. This seems to be the 
unlikely case. The first jet may generate a flash of 
lights when it breaks out the star (e.g., Ramirez- 
Ruiz, MacFadyen & Lazzati 2002); but the flash 
should have a thermal spectra, being inconsistent 
with the observed thermal spectra. 

Nevertheless our model, which is basically a 
standard jet-type model, can explain both the 
non-thermal spectra and the long delay as follows. 
As shown in Figure 1, the first jet is highly rela- 



tivistic and is traveling into the CSM. Such a jet 
can produces an external shock and non-thermal 
gamma-rays (e.g., Piran 2005). We propose that 
this is the observed precursor. 

The relatively soft spectra may be reproduced if 
T is not too large, although determining all other 
synchrotron parameters are not so simple. The 
small event count is consistent with this model, 
because the external shock colliding with smooth 
CSM is inefficient in emitting gamma rays. 

(vi) Now the last question is how the long delay 
is explained. The delay can be explained if the 
location of the external shock is separated from the 
closest internal shock by the distance of c x i^ei ~ 
3 x 10 12 cm in the local frame (Fig. 2 left-bottom). 
This condition appears to be difficult to satisfy 
because as mentioned in (iii) the GRB variability 
requires that the typical separation between two 
internal shocks is I ~ (3 - 60) x 10 10 cm. 

This apparent problem can be solved if we 
assume the presence of a sub-relativistic dense- 
shells somewhere behind the first shell. Such 
dense-shells may be formed from a shear-driven 
instabilities. Although the existence of relatively 
slow shells have not been studied much previously, 
some variations in the velocities of the shells are 
necessary for the internal shock model to work. 
This is because if all the shells are highly rela- 
tivistic, the kinetic to thermal energy conversion 
by the collisions of these shells will be quite inef- 
ficient (e.g., Piran 2005). 

Let us consider a relativistic shell with a mass 
M and the Lorentz factor T colliding with a rel- 
atively slow dense-shell. The relativistic shell will 
be decelerated to ~ T/2 when the shell collides 
with a mass of ~ M/(2r). Therefore, the shell 
is decelerated to sub-relativistic when the mass of 
the dense-shell is larger than ~ M. If the dense- 
shell has a mass much larger than M, relativistic 
shells are almost completely stopped, then large 
amount of thermal energy is generated and may 
turn into gamma-ray emission. 

If such a dense-shell is non-relativistic, then 
the distance between the shell and the " first" - 
relativistic shell increases with time (see Figure 
2 left). The main GRBs that occur behind the 
dense-shell will be visible only after the dense-shell 
becomes optically thin at t opt . The delay time is 
then estimated by idci — 'opt/c, where Z op t is the 
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distance between the dense-shell and the " first" - 
shell at t op t- 

For the model in Figure 1, at t = 100 sec, the 
"first" -shell is located at around radius = 94, in 
the units of c x 1 sec =1. Let the position, or the 
distance from the center of the explosion, of the 
shell at time t be x\, then it is approximately x\ ~ 
94 + c x (t — 100) = t — 6. Let us consider a dense- 
shell located at radius = .t 20 (^ x\ — I) at t = 100 
with radial velocity a(< 1). The position of the 
dense-shell at t is written as x 2 = x 2 o +a(t — 100). 
In terms of these quantities tdei is written as 
i do i ~ (xi - x 2 )/c = (1 - a)t opt - 6 + 100a - x 2a . 
Here we assume that the velocity of the dense- 
shell is constant, but it can increase with time by 
the interaction with surrounding material and the 
collision with the relativistic matter from behind. 
For simplicity let us assume that the shell has a 
constant velocity a until it becomes relativistic at 
time i re i. Then the above condition can be rewrit- 
ten as 

i d ci ~ (1 -a)t 2 -6 + 100a -:r 2 o, (2) 

where t 2 =min(t opt , t rc i)- F° r example, when 
x 20 — 80, a = 0.85, and tdei > 100, t 2 is con- 
strained as t 2 > 670 sec. 

There is another constraint that has to be sat- 
isfied: if the jet duration time, Atj ct , is too short, 
the collisions of shocks finish before the GRBs 
become visible. Let x% be the position of the 
"last" relativistic shell and suppose that x 3 — x 30 
at t = Aijet, then x 3 (t) ~ x 3a + (t - Ai jet ). 
Here typically x^o is x^q < cxlO sec. This con- 
straint is then written as x^(t 2 ) < x 2 (t 2 ), or, 
Aij et > ^30 - x 2Q + (1 - a)t 2 + 100a. Using equa- 
tion (2), this becomes Aij ct > x^q + t^ c \ + 6. If 
idol + £30 > 100, this constraint is roughly the 
same as Equation (1). 

We note that some GRBs do not have any indi- 
cations of precursors. In our model this is simply 
understood as the very small td e \- idci may be 
nearly zero if there is no slow dense-shell or if the 
dense-shells are accelerated to relativistic speed 
quickly by the collisions with following shells. 

4. Conclusion 

We find that the relativistic jets from a C+O 
star can reproduce observed GRBs with certain 



conditions. Two conditions are especially impor- 
tant. One is the sufficiently long duration time 
for the jet, Atj et > 100 sec. In this Letter we 
treat the central engine as a black box, thus being 
unable to answer if such a condition is actually 
satisfied. However, our model may explain such 
duration because HNe is likely to be an asphcri- 
cal explosion where large fall-back may occur from 
equatorial direction possibly for a long time. We 
will investigate such explosions in future works. 

The other important condition is the existence 
of a sufficiently massive sub- relativistic "dense- 
shell" . The initial mass and velocity of such a shell 
together with the sufficiently long Atj et can pro- 
duce long tdei, as well as short idei from a relatively 
short time scale i V ari — 1 — 20 sec. This conclu- 
sion is robust because the existence of relativistic 
shells in front of the slow-dense-shell merge into 
the first-shell without forming a strong internal 
shocks. In our model, if two or more slow dense- 
shells are formed, main GRBs arc split into two 
or more parts. This phenomena appear to be seen 
in some of actual GRB events. In order to answer 
how the dense-shells and other relativistic shells 
with mean separation / ~ 10 11 cm are formed, we 
at least need to perform higher resolution calcu- 
lations to study hydrodynamical instabilities. We 
also need to investigate the effects of He-envelopes 
for the formation and distribution of those shells. 
Other unconsidered effects such as magnetohydro- 
dynamical effects may be critical for the instabili- 
ties. We leave such studies for future works. 
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Fig. 1. Distribution of T along the jet axis for the model in Section 2, plotted against the radius, R, in the 
units of c x 1 sec =1 (Fig. la). The location of gamma-ray photosphere is shown by arrows. The relativistic 
shock becomes optically thin for t >^ 1600 sec. The energy contained in matter with T greater than the 
indicated values is shown in Figure lb. 
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Fig. 2. — Schematic pictures of our GRB model in the local frame (left), and gamma-ray flux in the observer's 
frame (right). The distance between the initially slow dense-shell and the front most shell, ' D' , determines 
the delay between the precursor and main-bursts. 
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